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Abstract. In earlier papers Saxena et al. (2002, 2003) derived the solutions 



> 

pij ' of a number of fractional kinetic equations in terms of generalized Mittag-Leffler 

\Q . functions which extended the work of Haubold and Mathai (2000). The object 

<^> of the present paper is to investigate the solution of a unified form of fractional 

kinetic equation in which the free term contains any integrable function f(t), 
which provides the unification and extension of the results given earlier recently 
^^ by Saxena et al. (2002, 2003). The solution has been developed in terms of the 

Q-i' Wright function in a closed form by the method of Laplace transform. Further 

f-H we derive a closed-form solution of a fractional diffusion equation. The asymp- 

totic expansion of the derived solution with respect to the space variable is also 
discussed. The results obtained are in a form suitable for numerical computa- 
tion. 



1 Introduction 

Fundamental laws of physics are written as equations for the time evolution of 
a quantity X(t), dX(t)/dt = —AXif), where this could be Maxwell's equations 
or Schroedinger's equation (if A is limited to linear operators), or it could be 
Newton's law of motion or Einstein's equations for geodesies (if A may also 
be a nonlinear operator). The mathematical solution (for linear operators) is 
X(t) = A(0)Exp{— At}. The initial value of the quantity at t = is given by 
X(0). 

The same exponential behavior referred to above arises if X(t) represents the 
scalar number density of species at time t that do not interact with each other. 
If one denotes A p the production rate and A^ the destruction rate, respectively, 



1 



the number density X(t) will obey an exponential equation where the coefficient 
A is equal to the difference, A p — Ad- Subsequently, A~ 1 is the average time 
between production and A^ is the average time between destruction. This 
type of behavior arises frequently in biology, chemistry, and physics (Kaplan 
and Glass, f998; Hilfer, 2000; Metzler and Klafter, 2000; Aslam Chaudhry and 
Zubair, 2002). This paper, in Section 3, considers the fractional generalization 
of the kinetic equation and derives closed form representation of its solution. 

The evolution of the number density X(x,t) relates the first derivative in 
time of this function to a spatial operator applied to the number density. The 
initial value of the function at time t — is given by X(x,0). In Section 4 of this 
paper, we give a relationship between the solution of the equation of evolution 
and the solution of the diffusion equation belonging to its fractional extension 
(Kaplan and Glass, f998; Hilfer, 2000; Metzler and Klafter, 2000). 

Section 2 summarizes mathematical results concerning solutions of the ki- 
netic and diffusion equations in Sections 3 and 4, respectively, widely distributed 
in the literature or of very recent origin. These involve the Mittag-Lefflcr func- 
tion, Wright function, and H-function, and the application of fractional calculus, 
Fourier transform, and Laplace transform to them. 

Specifically, Haubold and Mathai (1995) discussed a conjecture of a variation 
of the solar neutrino signal looking at solutions of the standard kinetic equation 
based on lifetime densities and Poisson arrivals. A closed form representation of 
the fractional kinetic equation and thermonuclear function is derived by Haubold 
and Mathai (2000) in terms of certain series which represents the Mittag- Lcfficr 
function (Mittag-Leffler, 1903, 1905). In order to extend the work of Haubold 
and Mathai (2000) three fractional kinetic equations are solved by Saxena et 
al. (2002) in terms of the generalized Mittag-Leffler functions in a closed form. 
Further extensions of these results are provided by a recent paper of Saxena et 
al. (2003). Section 3 and 4 of this paper provide a unification and extension 
of the aforementioned results on fractional kinetic equations by investigating a 
closed-form solution of a unified fractional kinetic equation in which the free 
term contains a function f(t) by the method of Laplace transform. The solution 
has been obtained in terms of the generalized Wright function in a closed form 
in Section 3. Also presented is a closed-form solution of a fractional diffusion 
equation in terms of an H-function. Its asymptotic expansion for large values of 
the space variable is discussed. The results derived in this paper are in a form 
suitable for numerical computation. 

2 Mathematical Prerequisites 

A generalization of the Mittag-Leffler function (Mittag-Leffler, 1903, 1905) 

£^)-E w Z " TV (aeC,i?e(a)>0) (1) 



was introduced by Wiman (1905) in the general form 

E «Az)--=J2 T (na + 0y («^eC,i?e(a)>0). (2) 

The main results of these functions are available in the handbook of Erdclyi, 
Magnus, Oberhettinger and Tricomi (1955, Section 18.1) and the monographs 
written by Dzherbashyan (1966, 1993). Prabhakar (1971) introduced a general- 
ization of (2) in the form 

*%,/»(*) = E tv {l) 2 Z L \v (o e ,P, 1 GC,Re(a)>0), (3) 

^ r (na + /?)(n)! 

where (7)0 = 1 is the Pochammer symbol, defined by 

( 7 )o = l,(7)fe = 7(7 + l)...(7 + fc-l) (fc=l,2,...), 7 ^0. (4) 

It is an entire function with p = [Re(a)]~ 1 (Prabhakar, 1971). 

For 7 = 1, this function coincides with (2), while for 7 = /? = 1 with (1): 

El (z) = E ai0 (z), E^z) = E a (z). (5) 

We also have 

#09,7,*)= iF 1 (p; r ,z) = T( 1 )E^(z) (6) 

where $(/3, 7; z) is Kummer's confluent hypcrgcomctric function defined in Erdclyi 
et al. (1953, p.248, eq.l). 

The Mellin-Barnes integral representation for this function follows from the 
integral 

F 7 M _i L f r(-pr(7 + 0(-zM m 

where u> = (— l) 1 / 2 . The contour fHs a straight line parallel to the imaginary 
axis at a distance 'c' from the origin and separating the poles of r(— £,) at 
the points £ = v [y = 0, 1, 2, . . .) from those of T(7 + £) at the points £ = 
—7 — 1/ (1/ = 0, 1, 2, . . .). If we calculate the residues at the poles of r(7 + £) at 
the points £ = —7 — v {y = 0, 1, 2, . . .), then it gives the analytic continuation 
formula of this function in the form 

vp ^ _ (~*)~ 7 f- r (^) (-*)"" 1 -1-1 m 

EaAz) "W h r[/J-a( 7 + „)] M ! ' ' z > L (8) 

From (8) it follows that for large z its behavior is given by 

ElJz)~0(\zD,\z\>l. (9) 



The H-function is defined by means of a Mcllin-Barnes type integral in the 
following manner (Mathai and Saxena, 1978, p. 2 ): 
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where 



©(0 = 



2ttw Jn 



e(0*-«de, (10) 



(11) 



Il?= m+ i r(i - 6, - BiOlin^i r( 0j + A-o] ' 

m,n 7 p,q G A^o with < n < p, 1 < m < q, Ai,Bj G R + ,ai,bj € i? or 
C(i = 1, . . . ,p; j = 1, . . . , q) such that 



Ai(bj + k) 76 Bj^oj - / - 1) (A;, / e iV ; i = 1, . 



,«;j 



(12) 



where we employ the usual notations: Nq — (0, 1,2,.. .); R = (—00, 00), R + = 
(0,oo) and C being the complex number field. Q is a suitable contour sepa- 
rating the poles of T(bj + sBj) from those of T(l — a,j — Ajs). A detailed and 
comprehensive account of the H-function is available from Mathai and Saxena 
(1978). 

It follows from (7) that the generalized Mittag -Leffler function ~EP a Az) can 
be represented in terms of the H-function in the form 



KA Z ) = 



1 



r(7) 



0-1,1 
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1(0,1), (l-fl.a) 



,(Re(a) > 0;a,/3,7 G C), 



and in terms of the Wright function as 
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(13) 
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where \^i(z) is a special case of Wright's generalized hypergcomctric function 
(Wright, 1935, 1940); also see Erdelyi et al. (1953, Section 4.1) where the 
function is defined by 
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where 1 + Y^j=i Bj ~ Y7j=i A? — 0' (equality only holds for appropriately 
bounded z). The relation connecting p ty q (z) and the H-function is given by 
Mathai and Saxena (1978, p. 11, eq. 1.7.8) 
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It is interesting to observe that for 7=1, (13) and (14) give rise to the following 
results for the generalized Mittag-Lefner function 



E a .A z ) 



where Re(a) > 0, a, (i, e C. 

If we further take f3 = 1 in (17) and (18) we find that 
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(18) 



£ a (z) = 1*1 

= H, 'n 



(1,1 
(1,«) 

■0 



(0,1) 
(0,l),(0,a) 



(19) 
(20) 
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where Re(a) >0,aeC. It is shown by Kilbas et al. (2003) that 

r- 1 (x - ^-^(u^^M* - t]")dt 

where p, fj,, j,a,u> € C, Re{v) > 0, Re(p) > 0, which is a generalization of the 
well-known result of Erdelyi et al.(1953, p. 271, cq. 6.10.15) 



t"- 1 (x - t)^ 1 <5(a, i>, w*)#(7, /U, u)[x - t])dt 
= B(v,fj,)x v+ ' i ~ 1 ®('y + cr,(J, + i';u!x), 



(22) 



where Re(p) > 0, Re{v) > 0., The Laplace transform of the H-function in terms 
of another H-function is given by Prudnikov ct al. (1989, p. 355, eq. 2.25.3) 
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where o > 0,Re(s) > 7 Re[p + <<7< m (^)] > 0, \argz\ < [7r/2]0,0 > 0; 

n p m q 

= E^- E A i + E s ^ E ^- 

j=l j=n+l .7=1 j=m+l 



From (23) it can be easily seen that 
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(23) 



(24) 



(25) 



where a > 0, i?e(s) > 0, i?e p + a?g< n { i^} > 0, |ar 5 z| < |tt0i, 0i > 
0, where 0i = — a, is defined in (24). 



From Prudnikov et al. (1989, p. 355, eq. 2.25.3.2) and Mathai and Saxena 
(1978, p. 49), it follows that the cosine transform of the H-function is given by 



tp- 1 cos(kt)H™> n 



p.q 
p\ rrn+l.m 
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(26) 



(Tr/n^i+a 



k> 1 
a 



(l-6„B g ),(i + £,f) 



J-p,^p,),\.2 TJ^Jj 



where Re[p + fx ^ n < m (^)] > 0, „ > 0, Re[p + M ^J^j 1 )] < 
1, \arg a\ < (?r0/2), 6 > 0, 9 is defined in (24). 

The Ricmann-Liouvillc fractional integral of order v G C is defined by Miller 
and Ross (1993, p. 45; see also Srivastava and Saxena, 2001) 



o^r/w = 



i 



(t - uY^fiujdu, 



(27) 
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where Re(v) > 0. Following Samko ct al. (1993, p. 37) we define the fractional 
derivative for a > in the form 



1 



d" 



f(u)du 



oD?f(t) = — / -— / , J \' — -r, (n=[Re(a)] + l), (28) 

* J w T(n-a) dt n J (t - u)*-^ 1 ' v L v yj h K ' 



where [i?e(a)] means the integral part of Re(a). 
In particular, if < a < 1, 



oD?f(t) 



1 



f(u)du 



dt T(l - a) J (t - w) Q 

and ii a = n E N — (1,2, . . .), then 

D?f(t) = D n f(t)(D = d/dt) 

is the usual derivative of order n. 

From Erdelyi et al. (1954, p. 182) we have 



L{ Drf{t)} = *- v m, 



where 



F(s)=L{f(t);s} 



l f{t)dt, 



(29) 
(30) 

(31) 
(32) 



where Re(s) > 0. 

The Laplace transform of the fractional derivative is given by Oldham and 

Spanier (1974, p. 134, eq. 8.1.3; see also Srivastava and Saxena, 2001): 



L[ D? f(t)} = s a F(s) - ]T s"- 1 Dr k f(t)\t= 



(33) 



fe=i 



In certain boundary-value problems arising in the theory of viscoelasticity and 
in the hereditary solid mechanics, Caputo (1969) introduced the following frac- 
tional derivative of order a > 0, defined by 

1 /"' f m (T)d,T 

D?f(t) = w 7/ , \ t\ (m-l<a<m),Re(a) >0,meiV 

* w r(m-a)J (t-T) a+1 - mK ~ ' v ' 

d m 
= _/(*), if a = m (34) 

Caputo (1969) has also shown that 

rn—l 

L{D?f(t);s} = S Q F( S ) - ^ s «- fc - 1 / fc (0+), (m - 1< a < m). (35) 

fe=0 

The formula (35) is very useful in the solution of diffcrintegral equations gov- 
erning certain physical problems. 

3 Unified Fractional Kinetic Equation 

If we integrate the standard kinetic equation, we obtain 

N,(t) - N = -c t Q D^ Ni(t), (Ci>0), (36) 

where oD^ is the standard Ricmann integral operator. In the paper of Haubold 
and Mathai (2000) the number density of the species i, Ni — Ni(t), is a function 
of time and Ni(t = 0) = Nq is the number density of species i at time t — 0. If 
the index i is dropped in (36), then the solution of the fractional kinetic equation 

N(t) - N = -c v D^N(t), (37) 

is obtained in a series form 

Applying the definition (1), the above equation can be rewritten in a compact 
form in terms of the Mittag-Leffler function 

N(t) = N E v [-(cty],v>0. (39) 

Theorem 1. If c > 0, v > 0, then for the solution of the integral equation 

N(t) - N f(t) = -c v Dt"N(t), (40) 

where f(t) is any integrable function on the finite interval [0, b] , there holds the 
formula 

N(t) = cNoj^H 1 ^ [c"{t - t) v |S:l;;j f(0 , ] f(.r)dr, (41) 



where H 1 : 2 (.) is the H-function defined by (10). 

Proof: Applying the Laplace transform to equation (40) and using (31) gives 



N*(s) = L[N(t);s] = N 

where F(s) is the Laplace transform of /(£); 
Since (Mathai and Saxena, 1978, p. 152) 



F(s) 
1 + (§)"' 



(42) 



~ - #1,1 (-) l(l.l) 



then using (25), we obtain 



L- 



H 



(7)"i8;l!]]=*- 1 ^ 1 [(ct)-'i|i;ij 



^1(1,1), (0,^)i 



(43) 



(44) 



By virtue of the following property of the H-function (Mathai and Saxena, 1978, 
p.4, eq. 1.2.2) 

(45) 
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= ciJ 



1.1 
1.2 



(ci) 1 



(-i.l) 
(-i,l),(0^) 



(46) 
(47) 



Eq. (47) follows from (46) if we use the formula (63). 

Taking the inverse Laplace transform of (42) by using (47) and applying the 

convolution theorem of the Laplace transform we arrive at the desired result 

(41). 

If we set f(t)~t p 1 , we obtain the following result established by Saxena 
et al. (2002, p.283, cq. 15) 

Corollary 1.1. If v > 0, p > 0, c > 0, then for the solution of the integral 
equation 

N(t) - Not"- 1 = -c v Dt V N(t), (48) 

there holds the formula 



N(t) = N T{p)tP- l E u , p [~{ctT] 



(49) 



Eq. (49) can be established by expressing the H-function, occurring in (41) in 
terms of an equivalent series by means of the computable representation of the 
H-function (Mathai and Saxena, 1978, p. 71, cq. 3.7.1.) and integrating term by 
term by means of the beta function formula. 

On the other hand, if we take f(t) = t p_1 -EJ„[— (ci)^], then another result 
recently given by Saxena et al. (2003) is obtained: 



Corollary 1.2. If c > 0, v > 0, p > 0, then for the solution of the integral 
equation 

N(t)-N t»- 1 E?J-(cty} = -c" Dr v N(t), (50) 

there holds the formula 

N(t) = N t»- 1 EJ+ 1 [-(cty}. (51) 

Eq. (51) can be proved with the help of the results (13) and (21). 

4 A Fractional Diffusion Equation 

In the following we derive the solution of the fractional diffusion equation using 
(52). The result is obtained in the form of the following 

Theorem 2. Consider the fractional diffusion equation (Metzler and Klafter, 
2000; Jorgenson and Lang, 2001) 

D»N(x, t) y^-J(x) = c^N(x, t), (52) 

with the initial condition 

D^ k N(x,t)\ t = = 0, (k = l,...,n), (53) 

where n = [i?e(z/)]+l, c v is a diffusion constant and 5(x) is Dirac's delta function. 
Then for the solution of (52) there exists the formula 



N(x,t) = — l - TT7 ^ H i2 



\A 2 ,(i-f^) 



Ac u t v 1(0,1), (1/2,1) 



(54) 



Proof. In order to derive the solution of (52), we introduce the Laplace- Fourier 
transform in the form 



oc />oo 



N*(k,s)= / e- st+lkx N(x,t)dxdt. (55) 

J0 J -co 

Applying the Fourier transform with respect to the space variable x and Laplace 
transform with respect to the time variable t and using (53), we find that 

s"N*(k, s) - s"- 1 = -c"k 2 N*(k, s). (56) 

Solving for N*(k, s) gives 

N*(k,s)= ^ ,„ . (57) 

y ' s v + c u k 2 K ' 



To invert equation (57), it is convenient to first invert the Laplace transform 
and then the Fourier transform . Inverting the Laplace transform , we obtain 



1/^2.1/1(0,1) 
C K l 1(0,1), (0,i/) 



N(k,t) = E v (-c v k 2 t v ), 
which can be expressed in terms of the H- function by using (18) as 

N(k,t) = H\* 

Using the integral 

1 r°° 1 f°° 

— / e- lkx f(k)dk = - f(k)cos(kx)dk, 

2tT J_ oc 7T J 

and (26) to invert the Fourier transform, we see that 



N(x,t) 



TV 



o 



cos(kx)H l ' 2 c k t \( \) (0,1/ 



dk 



(58) 



(59) 



(60) 



(61) 



\x\ 3 ' 3 



\X\ 2 ,(1,1),(1,„),(1,1) 
c v t v 1(1,2), (1,1), (1,1) 



Applying a result of Mathai and Saxena (1978, p. 4, eq. 1.2.1) the above expres- 
sion becomes 



N(x,t) 



1 



tH. 



2.2 



N_ i (1,1/), (1,1) 

c „ r 1(1,2), (1,1) 



If we employ the formula (Mathai and Saxena, 1978, p. 4, cq. 1.2.4): 



<T zj-m,n 



x°K 

Eq. (62) reduces to 

7V(x,t) 



[dp -Ap) 
(b q ,B q ) 



Tjm.n 



(a p +aA p ,A p ) 

(b q +aB q ,B q ) 



H: 



2.0 



(c^)l/2 2,2 



k| 2 ,(l-f ,i/),(l/2,l) 

c i/ r 1(0,2), (1/2,1) 



In view of the identity in Mathai and Saxena (1978, eq. 1.2.1), it yields 



1 



N fa *) = 7 777? H l'l 



c^i" 



(1-7.") 

(0,2) 



Using the definition of the H-function (10), it is seen that 



N(x,t) 



r(20 



27r W (c"i-)i/ 2 A, T[l - | + */£] 



c^i 7 - 



-s c 



«. 



(62) 



(63) 



(64) 



(65) 



(66) 
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Applying the well-known duplication formula for the gamma function and in- 
terpreting the result thus obtained in terms of the H-function, we obtain the 
solution as 



N(x,t) 



1 



--H 



x 



(0,1), (1/2,1) 



(67) 



Finally the application of the result (Mathai and Saxena (1978, p. 10, eq. 1.6.3) 
gives the asymptotic estimate 



N(x,t)~0{ [M^][exp. 



(2-^)(H 2 Q 



{0<is< 2). (68) 
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